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■ ABSTRACT 

. We present long-slit observations in the optical and near infrared of eight H n 

regions in the spiral galaxy NGC 4258. Six of the observed regions are located in the 
I , SE inner spiral arms and the other two are isolated in the northern outer arms. 

■ A detailed analysis of the physical conditions of the gas has been performed. For 
(yT) ' two of the regions an electron temperature has been derived from the [SIII] A 6312 A 
Q\ , line. For the rest, an empirical calibration based on the red and near infrared sulphur 
t— I ■ lines has been used. The oxygen abundances derived by both methods are found to 

be significantly lower (by a factor of two) than previously derived by using empirical 
calibrations based on the optical oxygen lines. 

In the brightest region, 74C, the observation of a prominent feature due to Wolf- 
Rayet stars provides an excellent constraint over some properties of the ionizing clus- 
• ters. In the light of the current evolutionary synthesis models, no consistent solution is 

Q-f found to explain at the same time both the WR feature characteristics and the emis- 

1 , sion line spectrum of this region. In principle, the presence of WR stars could lead to 

■ large temperature fluctuations and also to a hardening of the ionizing radiation. None 
' of these effects are found in region 74C for which the electron temperatures found from 

the [SIII] A 6312 A line and the Paschen discontinuity at 8200 A are equal within the 
' errors and the effective temperature of the ionizing radiation is estimated at around 

; 35300 K. 

. Both more observations of confirmed high metallicity regions and a finer metallic- 

' ity grid for the evolutionary synthesis models are needed in order to understand the 

ionizing populations of HII regions. 

Key words: galaxies: NGC 4258 - galaxies: HII regions - galaxies: stellar content - 
stars: WR stars 



1 INTRODUCTION 

Theoretically, the evolution of a young stellar population 
depends on metallicity at least through two very important 
effects: the increasing opacity of the stellar material and the 
dependence of mass loss on metal content in high mass stars. 
As a consequence of the first, the effective temperature of 
ionizing stars should be lower in regions of higher metallicity 
(see for example McGaugh 1991) but, on the other hand, as 
a consequence of the second, that is if the strength of stellar 
winds increases with metallicity, the loss of the outer en- 
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velopes of the most massive stars can increase their surface 
temperature to very high values. These highly evolved mas- 
sive O stars are identified with the Wolf Rayet population. 

Most theoretical evolutionary models for ionizing star 
clusters predict the appearance of WR stars - although at 
slightly different cluster ages depending on the specific de- 
tails of the models - and in all of them the fraction of WR 
stars increases with metallicity since the limiting mass for a 
star to enter the WR phase decreases with increasing metal- 
licity. However, the effect of these stars on their surround- 
ing ionized gas depends on the characteristics of the as- 
sumed wind opacity and therefore while some models predict 
the existence of high metallicity Hll regions of high excita- 
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tion (e.g. Garcia- Vargas & Diaz 1994; Stasiriska & Leitherer 
1996), others do not (e.g. Cervino & Mas-Hesse 1994). The 
finding of this extreme HII region population would consti- 
tute a firm evidence of the existence of WR stars with high 
effective temperatures, thus favoring the first kind of models 
but, even if these regions are not definitively identified, the 
study of high metallicity Hll regions is of crucial importance 
to investigate the processes of star formation and evolution 
in high metallicity environments. This investigation gains 
even more relevance if we take into account the recent ob- 
servations that point out to a firm connection between star 
formation and activity in active galactic nuclei (Heckman 
et al. 1997, Gonzalez-Delgado et al. 1998) already predicted 
on theoretical grounds (Terlevich & Melnick 1985; Terlevich 
et al. 1992), since these nuclei seem to reach rather high 
metallicities (Phillips et al. 1984; Gonzalez-Delgado & Perez 
1996) as deduced from the study of their circumnuclear star 
formin g regions. 



Un fortunately, the term "high metallicity Hll region" is 
rather ill defined. Most authors apply the term to Hll re- 
gions of solar or oversolar metal content, but the difficulty 
to derive directly the metallicity of these cool HII regions 
is well known and the empirical calibrations commonly em- 
ployed provide estimations with an uncertainty larger than 
a factor of 2 (see e.g. Diaz et al. 1991). On the other hand, 
the characteristics of WR features are highly dependent on 
metallicity. Both the luminosity of the WR "bump" , that is 
the combination of the NlII and the Hell lines at around A 
4660 A, and its equivalent width are predicted to be largest 
at the highest metallicity and almost an order of magni- 
tude difference is found between the computed WR "bump" 
intensity relative to H/3 in clusters of half solar and solar 
metallicities (Schaerer & Vacca 1998). Therefore, the detec- 
tion of WR features combined with a detailed analysis of 
the emission line spectra of suspected high metallicity Hll 
regions, can constitute excellent tools for the simultaneous 
determination of both the metallicity and the age rather 
accurately in these regions. 

In the frame of a long term programme to study high 
metallicity Hll regions we are presently carrying out a search 
for the high metallicity and high excitation Hll region pop- 
ulation by performing spectrophotometric observations at a 
resolution high enough to detect and measure WR features. 
For that, we have selected from the literature the giant HII 
regions with solar or oversolar abundance, as deduced from 
empirical calibrations based on the optical forbidden lines, 
showing the highest excitation. This corresponds to a ratio 
of [OIII]/H3 ~ 1.0, which is actually quite moderate but 
still high compared to that expected in high metallicity HII 
regions. Here we present the first of these investigations in- 
volving several Hll regions in the galaxy NGC 4258. 

NGC 4258 (M106) is c lassif ied as SAB(s)bc and Sb(s)II 
t ypes by De Vaucouleurs (1976) and Sandage & Tammann 



(1981) respectively. Heckman (1980) and Stauffer (198S 



classified its nucleus as belonging to a transition type be- 
tween LINER and H II-region like emiss ion-line galaxy. As 
it is addressed by Courtes et al. (1993) (hereinafter C93), 



the H II region distribution is very peculiar in the sense 
that three general structures are found: the normal spiral 
arms, with dense and very bright Hll regions following pure 
spiral shapes as a consequence of the barred spiral type of 
this galaxy; the faint outer arms (H II region poor), and the 



anoma lous spiral arms, discovered by Courtes & Cruvellier 
(1961), wit hout star formation and w ith no evidence of very 
blue stars (Deharveng & Pellet 1970). It is assumed that the 



anomalous arms were formed by powerful jets ejected by the 
active nucleus of the galaxy. Stiiwe, Schulz & Hiihnermann 
showed that the active nucleus must be obscured 



( [1992D 



along the line of sight as is also indicated by its i nvisibil- 
ity in IUE spectra. Cecil, Morse & Veilleux (1995) found, 



from the emission-line flux ratios, LINER-like gaseous ex- 
citation along the jets. They also point out the possibility 
that the circumnuclear region has undergone recent forma- 
tion of massive stars induced by the jets that are flowing 
into the disk. However, no evidence for this star formation 
has been found, maybe due to the obscuration by molecu- 
lar gas along the line of sight. Finally, a subparsec diameter 
gaseous disk shows keplerian motion around the centre of 
the galaxy, which indicates the presence of a mass of 3.6 x 



10 M q within a radius of less than 0.13 pc (Miyoshi et al 



1995). This is probably the best existing evidence for the 
presence of a massive black hole at the centre of a galaxy. 

The problem of the distance determination for this 
galaxy still remains unsolved. The quoted distanc e values 
range from 3.4 Mpc (Roy, Arscnault fe Joncas 198(: ) to 10.4 
Mpc ( Sandage fc Tammann 1981|). The adop ted value in the 
present work is 5.5 Mpc (Martin et al. 1989). 



We have observed eight Hll regions along two different 
spiral arms. Four of them have been identified as: 74C, 69C, 
5N(A) and 5N(B) in C93. Regions 5N(A) and 5N(B) are lo- 
cated in the northern outer arm at a galactocentric distance 
of 20.3 Kpc. The other six regions are located in the SE in- 
ner arm at 8.05 Kpc from the gala xy centre. Region 74 C is 
a discrete radiosource at 4.9 GHz ( Hummel et al. 1989] ). 

Regions 74C a nd 5N were previously observed by Oey 
& Kennicut (|l993[) who derived values of 12+log(0/H) of 



8.87 and 8.74 for each of them respectively (the solar value 
is 8.92). With a ratio of [Olll] to H/3 of 0.98 and 1.44 both 
regions fit our selection criteria. 



2 OBSERVATIONS AND REDUCTIONS 

Our spectrophotometric observations were obtained with 
the 4.2m William Herschel Telescope at the Roque de los 
Muchachos Observatory, in 1994 March 13, using the ISIS 
double spectrograph, with the TEK1 and a EEV CCD de- 
tectors in the blue and red arm respectively. Two gratings 
were used, R300B in the blue and R316R in the red arm, 
covering from [O II] A3727 to [S III] A9532 in three different 
spectral ranges of ~1700 A each. The dispersion of 1.4 A 
pixel -1 with a slit width of l'/5 gives a spectral resolution 
of ~4 A. The nominal spatial sampling is 0'.'7 pixel -1 . The 
seeing was ~l'.'5. The slit was centered at the positions of 
regions 74C and 5N, as given by Oey & Kennicutt (1993). 
A third position was observed in order to include other HII 
regions outlining the SE inner arm. In all cases the position 
angle was chosen as to maximize the number of HII rgions 
in the slit. A journal of observations is given in Table 1. 

The data were reduced using the IRAF (Image Reduc- 
tion and Analysis Facility) package following standard meth- 
ods. The two-dimensional wavelength calibration was accu- 
rate to 1 A in all cases. The two-dimensional frames were flux 
calibrated using two spectroscopic standard stars observed 
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Table 1. 


Journal of observations 




P.A. 




Grating 


A range (A) 


Exposure (s 


255 


R300B 


3675-5300 


1800 


255 




R300B 


3675-5300 


1800 


255 




R316R 


5890-7600 


1200 


255 




R316R 


5890-7600 


1800 


255 




R316R 


7960-9670 


1800 


255 




R316R 


7960-9670 


1800 


150 




R300B 


3675-5300 


1800 


150 




R300B 


3675-5300 


1800 


150 




R300B 


3675-5300 


1200 


150 




R316R 


5890-7600 


1800 


150 




R316R 


7960-9670 


1800 


150 




R316R 


7960-9670 


1800 


133 




R300B 


3675-5300 


1800 


133 




R300B 


3675-5300 


1200 


133 




R316R 


5890-7600 


1800 


133 




R316R 


7960-9670 


1800 



a All the observations made on 1994 March 13/14 

during the same night with a 5'.' width slit. The agreement 
between the individual calibration curves was better than 
8% in all cases. The spectra were previously corrected for 
atmospheric extinction using a mean extinction curve appli- 
cable to La Palma observing site. The removal of the atmo- 
s pheric water-vapour absorp tion bands in the near infrared 
(Diaz, Pagel & Wilson 1985), was achieved dividing by the 
relatively featureless continuum of a subdwarf star observed 
in the same night as the galaxy. 

The most critical step in the reduction process was the 
background subtraction in the near infrared frames. There 
are numerous OH emission lines of the night-sky spectrum 
that contaminate the near infrared spectra. This contami- 
nation is particularly important in the range from 9376 to 
9600 A flOsterbrock, Fulbright fc Bida 1997| ). Hence, the 
[Sill] A9532 A and Pa 8 lines intensities could be affected 
even after the background subtraction is made. Fortunately, 
the [Sill] A9069 A and Pa 9 lines are unaffected by these 
night-sky lines. The observed A9532/A9069 ratio is above 
the theoretical value of 2.44 in four of the observed regions. 
In these cases, the [Sill] A9532 line was scaled by means of 
the more accurate [Sill] A9069 line intensity measurement. 



3 RESULTS 

Figure 1 shows the spatial distribution of the Ha flux along 
the slit for the three different positions observed. Regions 
5N A and B are clearly identified on the Ha profile cor- 
responding to PA=255°. Also regions 74C and 69C, show- 
ing a certain degree of substructure, are resolved at posi- 
tion angle PA=133°. Four different regions are resolved at 
PA=150°that we have named GA1, GA2, GA3 and GA4. 
A fifth region is clearly seen in the Ha profile that is not 
detected in the near IR frames. 

The spectra corresponding to each of the identified re- 
gions is shown in Figs. 2 to 4. WR features around A 4680 
A are seen in the spectrum of region 74C and to a lesser 
extent in that of region 5N. 



5N(A) 



5N(B) 




(b) P.A. 133" 



1.5X10- 16 - 



£. 10-'" 




Figure 1. Ho profiles for the three observed slit positions includ- 
ing: (a) regions 5N(A) and 5N(B); (b) regions 74C and 69C; (c) 
regions GA1, GA2, GA3 and GA4. 
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Figure 2. Merged spectra for the regions 5N(A) and 5N(B) 



3.1 Line intensities 

Emission line fluxes were measured using the IRAF SPLOT 
software package, by integrating the line intensity over 
a local fitted continuum. The errors in the line fluxes 



have been calculated from the expression oi = cr c N 1 ^ 2 [l 
+ EW/(NA)] 1/2 , where g i is the error in the line flux, a c 
represents the standard deviation in a box near the mea- 
sured emission line and stands for the error in the continuum 
placement, N is the number of pixels used in the measure- 
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Figure 3. Merged spectra for the regions GA1 and 69C. They could correspond to different regions of the same nebula 



merit of the line flux, EW is the line equivalent width, and 
A is the wavelength dispersion in angstroms per pixel. 

The observed line intensities relative to the H/3 line 
were corrected for interst ellar reddening a ccording to an 
average extinction curve jOsterbrock 198S ) and assuming 



the Balmer line theoretical values for case B recombination 



( Brocklehurst 1971). The presence of an underlying stellar 
population is evident in the blue spectra of regions GA2, 
GA3 and GA4. The H7 and H<5 Balmer lines are clearly af- 
fected by this stellar absorption. An iterative process was 
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Figure 4. Merged spectra for the regions GA2 and GA4 



applied in order to fit observed and theoretical Balmer line 
intensities and to obtain the reddening constant c(H/3). In 
all cases, Ha, H7 and H<5 were fitted, except in regions GA3 
and GA4 for which H<5 seems to be underestimated. Redden- 



ing corrected Paschen lines, when measured, are consistent 
with their theoretical values. 

The reddening corrected line intensities together with 
their corresponding errors are given in Table 2 for regions 
74C and 69C and Table 3 for the rest of the regions. Also 
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Figure 5. Merged spectra for region GA3. The bottom panel is an enlargement to show the prominent absorption features. 



given in the tables are the extinction corrected Ha flux and 
the H/3 equivalent width. For all the regions, Ha luminosi- 
ties are lower than 10 38 erg s -1 , with the exception of region 
74C, in whic h a value of 1Q 39 - 3 is found. According to Ken- 
nicutt (1983), this region can be classified as a supergiant 
Hll region. 



3.2 Physical conditions of the gas 

Electron densities for each observed region were derived from 
the ratio of the [Sll] AA 6717, 6731 lines following standard 
methods (see for example Osterbrock 1989). For regions 
GA2 and 74C densities of 160 cm -3 and 60 cm -3 respec- 
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Figure 6. Merged spectra with two intensity scales for region 74C. Notice the WR feature at 4686 A 



tively are obtained. For the rest of the regions the ratio of 
the [Sll] lines implies n e < 40 cm -3 . 

For two of the observed regions, 74C and 69C, it has 
been possible to measure the weak auroral [Sill] A 6312 A 
line, which, together with the observed intensities of the neb- 
ular lines at AA 9069, 9532 A, can be used to obtain a value 



for the electron temperature. Using the expressions given 
by Osterbrock (1989) and the effective collision strengths 
recently calculated by Tayal (1997) we find for t(S 2+ ) val- 
ues of 7500 ± 300 K and 7600 ± 500 K for regions 74C and 
69C respectively. 

Photoionization models indicate that, for abundances 
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Table 2. Reddening corrected line intensities for P.A.= 133° 



Region 




7 AH 
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( yo 


TT 1 n 
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HA _l_ 1 
Z<1 ± 1 
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A OOO 
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c(H/3) 




0.48 


0.08 


P(Ha) [J 




5576 


321 


EW(H/3)(A) 




81 


113 



a Affected by a cosmic ray 

6 10~ 16 erg cm -2 s _1 , corrected for reddening 



close to solar, the region where sulphur is twice ionized 
overlaps both the regions of once and twice ionized oxy- 
gen and therefore t(S 2+ ) is expected to be intermediate be- 
tween t(0 + ) and t(0 2+ ). Garnett (1992) from single star 
photoionization models finds a linear relation between the 
temperatures of 2+ and S 2+ 



t(0 



2 +1 



1.20i(S 2 



0.20 



(1) 



We have used that relation to obtain t(0 2+ ). Finally, 
Stasihska's (1980) relation between t(0 2+ ) and t(0 + ) has 
been used to calculate the temperature of the t(0 + ) zone. 
These different temperatures for regions 74C and 69C are 
listed in Table 4. 

Another independent way of temperature determina- 



tion is through the measure of the Paschen discontinuity at 
A 8200 A. The ratio of the Paschen discontinuity to the H/3 
intensity is, after correcting for the derived extinction value 
of c(H/3) = 0.48, 



APa 
I(H/3) 



(6.9 ± 0.8) x 10~ 15 Hz~ 



(2) 



Acc ordin g to the models computed by Gonzalez-Delgado et 
al. (1994), for the helium abundance derived for the region 



74C (see next section), this A Pa/I(H/3) implies an electron 
temperature of 8000 ± 700 K, which agrees with that de- 
rived from the [Sill] lines within the errors. Therefore, no 
temperature fluctuations are apparent in the region which 
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Table 3. Reddening corrected line fluxes for P.A.= 255°and 150° 
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might have been expected from the presence of WR stars 
(Perez 1997). 

For the rest of the regions it was not possible to obtain 
a direct measure of the electron temperature. An average 
temperature has been adopted from an empirical calibra- 
tion by means of the sulphur abundance parameter S23 (S23 
= ([Sn]6717, 6731 + [Sm]9069, 9532)/H/3) (Diaz & Perez- 
Montero 1999). The values for these temperatures, together 
with similar ones derived for regions 74C and 69C, are also 
given in Table 4. 



3.3 Chemical abundances 

Ionic abundances of oxygen, nitrogen and sulphur have been 
derived following standard methods (Pagel et al. 1992) and 



using the temperatures found above. These abundances are 
listed in Table 4. We have assumed that most of the oxygen 
and sulphur are in the first and second ionization stages and 
therefore O/H = 0+/H+ + 0++/H+ and S/H = S+/H+ + 
S ++ /H + . This assumption seems to be justified given the 
relatively low estimates of the electron temperature. 

The helium abundance was determined from the Hel A 
4471 and 6678 A lines. Mean values of 0.070 ± 0.005 and 
0.060 ± 0.005 respectively are obtained for regions 74C and 
69C. We have estimated the contribution of neutral helium 
from the expression: 



He° + 



He+ 



1 - 0.25 



0+ 
O 



He+ 
H+ 



(3) 



(Kunth & Sargent 1983). The correction factors obtained for 
regions 74C and 69C are 1.17 and 1.20 respectively which 



yield He abundances y=0.081 for region 74C and y=0.068 
for region 69C. 



3.4 Wolf-Rayet features 

A prominent Wolf-Rayet feature has been observed in re- 
gion 74C (se e Figure 6). Assum ing a distance to NGC 4258 
of 5.5 Mpc (Martin et al. 1989) and a constant extinction 
value through this region of c(H/3) = 0.48, as derived from 
the Balmer and Paschen recombination lines, the total lumi- 
nosity of the broad Hell feature, without the contribution of 



the [Fein] A 4658 A line, is (1.1 ± 0.1) x 10 38 erg s" 1 . This 
value comprises the broad features of both NlII AA 4634, 
4640 and Hell A 4686 A lines. The contribution of the NlII 
lines to the W R bump is metallicity-dependent according to 
Smith (1991). In our case this contribution represents 0.4 



times the total emission; hence we obtain for the Hell line 
luminosity: 



L(#e77A4686) 



(6.6 ± 0.6) x 10 37 erg s" 1 



(4) 



The non-detection of Nv AA 4604,4620 A, suggests that 
the WR stars are of late N or interm ediate type. Using the 
calibration of Vacca8& Conti ( |1992[ ), 38 ± 5 WN stars are 
found in this region. 

Another fainter Wolf-Rayet feature has been observed 
in region 5NA. The bump luminosity is 2.2 x 10 3B erg s -1 . 
The NlII lines contribution represent 0.45 times the total 
emission, hence the Hell line luminosity is 1.2 x 10 36 erg 
s -1 , which is compatible with the presence of 1 WN star. 
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Table 4. Derived physical conditions of the gas in the observed HII regions 
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Figure 7. Wolf-Rayet feature at A 4686 A in region 74C 

4 FUNCTIONAL PARAMETERS OF THE 
OBSERVED HII REGIONS 

Three are the fundamental paramet ers which contro l the 
emission line spectra of Hll regions (Diaz et al. 1991): the 
ionization parameter, the shape of the ionizing continuum, 
and the metallicity. 

The ionization parameter - i. e. the ratio of the ionizing 
photon density to the particle density - is a measure of the 
degree of ionization of the nebula and can be deduced from 
the ratio of two lines of the same element corresponding to 
two different ionization states, e.g. [Oll]/[OlIl] or [Sll] / [Sill] . 
Alternatively, it can also be determined from [Oll/H/3] or 
[Sll/Ha] if the metallicity of the region is known (Di'az 1994). 
By using a large grid of single star photoionization models 
the following expressions are found (Diaz et al. 1991, Diaz 
1999): 

logU = -l.mog{[OII]/Hf3) + 0.87log(Z/Z e ) - 1.68 (5) 
logU = -1.4Qlog([SII]/H/3) + l.Mog(Z/Z @ ) - 3.26 (6) 
logU = ~0miog{[On]/[OIII]) - 3.02 (7) 
logU = ~imiog([SII]/[SIII}) - 2.99 (8) 



We have derived U from the four expressions above. In 
all cases the value of U derived from the [Oll]/[OlIl] ratio is 
systematically lower than the rest thus implying low effec- 
tive temperatures for the ionizing stars (see e.g. Diaz 1999). 
We have therefore discarded this value and computed U as 
the mean of the other three. These adopted ionization pa- 
rameters - logU - are listed in Table 5 and their uncertainty 
is estimated to be around ± 0.2 dex. 

The shape of the ionizing continuum is directly related 
to the effective temperature of the stars that dominate the 
radiation field responsible for the ionization of the nebula. A 



recent ver sion of the photoionization code CLOUDY (For- 



land 1996j ) has been used to estimate the mean effective 



temperature of these stars. We have used Mihalas NLTE 
single-star stellar atmosphere models, with a closed geome- 
try and a constant particle density through the nebula. For 
all the regions, consistency is found for effective tempera- 
tures around 35,000 K - 36,000 K. The results from these 
photoionization models compared to the observations are 
presented in Table 5. 



5 DISCUSSION 

The Ha fluxes for regions 74C and 69C, uncorrected for 
reddening, are 2692 and 287 x 10 -16 erg s _1 cm -2 respec- 



tively. These fluxes are to be compared with 3492 and 332 x 
10~ 16 erg s _1 cm -2 measured by C93. Considering that our 
observations have been obtained through a narrow slit and 
region 74C is rather extended, the agreement can be consid- 
ered satisfactory. The combined Ha flux of our regions 5NA 
and 5NB is 187 x 10~ 16 erg s -1 cm -2 close to the value 
of 249 x 10~ 16 erg s" 1 cm' 2 given by C93 for region 5N. 
Regarding the rest of the regions, we tentatively identify re- 
gions GA2, GA3 and GA4 with regions 59C, 58C, 54C of 
C93. Region GA1 probably corresponds to the outer parts 
of region 69C. 

For each of our observed regions we have calculated the 
Ha luminosity and the number of hydrogen ionizing pho- 
tons from the observed Ha flux corrected for their derived 
reddening (see Tables 2 and 3). Both quantities depend on 
the distance D to NGC 4258 according to the expressions: 



L(Ha) = 3.62 x 10 3 



F(Ha)\ ( D 



erg s 



Q(H) = 2.65 x 10 4 



/F(Hc 



10- 



D_ 



photons s 



(9) 



(10) 
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where F(Ha) is expressed in erg s 1 cm 2 and D is given in 
Mpc. 

These values are given in Table 6 for the adopted dis- 
tance of 5.5 Kpc. Only region 74C has an Ha luminosity 
greater than 10 39 erg s _1 and can be classified as a super- 
giant HII region as defined by Kennicutt (f983). The rest 
of the regions have Ha luminosities typical of HII regions in 
early spiral galaxies, although all of them are greater than 
10 37 erg s _1 (Q(H) > 10 49 photons s _1 ) requiring more than 
a single star for their ionization (Panagia 1973). 

Filling factors for each observed region can be deter- 
mined from the reddening corrected Ha flux, F(Ha), and 
the derived ionization parameter, U, according to the ex- 
pression: 



0-34 (|P$) 



-1/2 



(5.5) (io- a ) 



3/2 



(11) 



( c B (H°,T) \ 1 I n c \ 
\y 10~ la J \W0) 



-1/2 



where as(ii"°,T) is the recombination coefficient for hydro- 
gen and n e is the electron density of the emitting gas. The 
fillinO factors computed in this way are listed in Table 6. We 
have used a value of a B (H°,T) = 3.76 x ICT 13 cm 3 s~\ 
corresponding to T= 7000 K and n e — 100 cm" 3 (Osterbrock 
1989). For regions 74C and GA2 we have used values of n e = 
60 and 160 cm -3 respectively, as derived from the ratio be- 
tween the red [SII] lines. For the rest of the regions for which 
only upper limits to the density could be derived, a value of 
n e = 10 cm -3 has been assumed. 

According to C93 the angular sizes of regions 74C, 69C 
and 5N are 38.6, 16.8 and 20.5 arcsec respectively. The corre- 
sponding linear sizes, at the assumed distance to NGC 4258, 
are 1.0, 0.45 and 0.55 Kpc . These are comparable to the 
sizes of Hll regions of moderate intensity found in other spi- 
ral galaxies (e.g. Gonzalez-Delgado et al. 1997). 

Ha angular effective diameters - that is the diameters 
containing half the Ha emission - of regions 74C, 69C and 
5N are 4.3, 3.7 and 6.1 arcsec respectively (C93) which trans- 
lates into 115, 99 and 163 pc. From the definition of U it is 
possible to obtain the angular size of the emitting region us- 
ing the reddening corrected Ha flux and the derived electron 
density: 



10" 



= 0.06 



( F ( Hc 



10- 



1/2 



u 

io- : 



-1/2 



(—) 

\ 100/ 



-1/2 



(12) 



where (j> is the angular diameter in units of 10 arcsec . This 
angular diameter does not depend on the assumed distance 
to the galaxy. For region 74C the derived angular diameter 
is 3.3 " close to the value of 4.3 " given by C93. For region 
69C their observed value of the effective diameter of 3.7 " 
yields a value of the electron density of n e = 5 cm -3 close 
to the assumed one of 10 cm -3 . Using this value of n e for 
the rest of the observed regions, except GA2 for which n e 
= 160 cm' 3 we obtain angular diameters between 1 and 2.6 
arcsec. 

We can also calculate the mass of ionized hydrogen fol- 
lowing the expression: 



M(HII) = 



m p Q(H) 



(13) 



n e (l + y+)a B (H°,T) 
(Osterbrock 1989) where m p is the mass of the proton and 
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Figure 8. Spectral energy distribution of different ionizing clus- 
ters for region 74C 



y the abundance by number of ionized helium, that we have 
taken as 0.10. Using observed and derived quantities this 
expression becomes 



M (m) . 53 s(^))(i)-(iL)V, 



(14) 



The corresponding masses of ionized hydrogen in each ob- 
served region are given in Table 6. 

From evolutionary models of single ionizing clusters and 
radiation bounded Hll regions, a relation between the num- 
ber of ionizing Lya photo ns per seco nd per solar mass and 
the H/3 equivalent width ( Diaz 199s| ) can be found, which 
allows the estimation of the ionizing cluster mass by taking 
into account the cluster evolution. 



log{Q{H)/M Q ) = 0.86log(EW(H(3)) + 44.48 



(15) 



For our observed regions all values for Q(H) are around 10 49 
photon s _1 except for region 74C for which a value of the 
order of 10 51 photon s _1 is found. Hence, in the absence 
of dust, a lower limit for the mass of the ionizing clusters 
can be estimated by means of the H/3 measured equivalent 
width and the Ha luminosity for each region. The estimated 
ionizing cluster masses range from 4 x 10 3 Mq for region 
GA4 to 1.12 x 10 5 M for region 74C. 

The high signal-to-noise spectrum of this latter region 
allows a more detailed modelling. The presence of WR fea- 
tures, in principle, provides a means to constrain the age 
of the ionizing population. The observed WR "bump" lu- 
minosity relative to H/3 (L(WR)/H/3) is 0.15. According to 
the models by Schaerer & Vacca (1998), at metallicity Z= 
0.008 (0.4 solar), these high values are only found at an 
age of 4.5 Myr. However, this ratio is strongly dependent 
on metallicity and changes by more than an order of mag- 
nitude between Z=0.008 (0.4 solar) and Z= 0.02 (solar). 
Therefore, the slightly higher metallicity of region 74C (Z 
~ 0.012) might allow earlier ages for the WR stellar popu- 
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Table 6. Physical properties of the observed HII regions 
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lation. Clearly, the available grid for stellar evolution mod- 
els is too coarse to elucidate this matter. The same applies 
to the equivalent width of the "bump" (EW(WR)) whose 
observed value is 10.8 A. Since the NIII A 4640 A line is 
strongly dependent on metallicity (Smith et al. 1996) it is 
probably more advisable to use only the flux and equivalent 
width of the Hell line. In that case, our observed values of 
L(HeII)/H/3 and EW(Hell) are 0.09 and 6.6 A respectively. 
Again, consistency is found for a model of Z=0.008 and an 
age of 4.5 Myr. 

Since, according to the recent evolutionary models of 
Leitherer et al. (1999; STARBURST99) , the observed value 
of the equivalent width of H/3 (81 A) points also to an age of 
around 4.5 Myr, we have computed the emission line spec- 
trum corresponding to this single ionizing cluster. Q This 
set of models makes use of the same evolutionary tracks as 
the models of Schaerer & Vacca (1998) and therefore pro- 
vides a selfconsistent frame. The results of the computation 
can be seen in Table 7 together with the observed emission 
line ratios. It can be seen that the model ionizing spectrum 
results too hard and does not reproduce the observations. 
This is not surprising in view of the low effective temper- 
ature found from single star photoionisation models (Table 
5). It can be better appreciated in Fig 8 where the spectral 
energy distribution of the cluster can be compared to that of 
the Mihalas stellar atmosphere of a star with T e / / = 35300 
K which adequately reproduces the data. A high energy tail 
is clearly present in the ionising cluster as the result of the 
WR star contribution. 

These models assume an enhanced mass loss during the 
evolution of massive stars, as prescribed by the stellar evo- 
lution models of Meynet et al. (1994). More moderate mass 
loss rates are assumed in the models of Bressan et al. (1993) 
and Fagotto et al. (1994) and with them a single population 
of 5.3 Myr can adequately reproduce the emission line spec- 
trum (see Table 7). The corresponding ionizing spectrum is 
also shown in Fig 8. However, this ionizing cluster provides 
only a few WR stars (about 3x 10 -5 per solar mass of the 
ionizing cluster, for a metallicity Z=0.008; Garcia- Vargas, 
Bressan & Diaz 1995). This is probably not enough to re- 
produce the observed values of L(WR)/H/3. 

The WR features observed in region 74C are similar 
to those found in other extragalactic HII regions. Table 8 
shows the values of the WR feature intensities and equiva- 
lent widths for 74C and other three well studied extragalac- 

t All the models are computed assuming a standard Salpeter 
IMF for the stellar ionizing cluster and using the photoionization 
code CLOUDY (Ferland 1996). 



tic HII regions: NGC 604 in M33 (Diaz et al 1987; Terlevich 
et al. 1996), region A in NGC 3310 (Pastoriza et al. 1993) 
and region A in NGC 7714 (Garcia- Vargas et al. 1997). Also 
given in the table are the metallicities of the regions in the 
form 12+log(0/H) and the effective temperature of the ion- 
izing radiation estimated from single star photo-ionization 
models. 

The tabulated values are difficult to explain in the 
light of the recent models of Schaerer & Vacca (1998). 
The expected decrease of the intensities of WR features 
with decreasing metallicity is not clearly observed. How- 
ever, a decrease of the temperature of the ionizing radia- 
tion with increasing metallicity is apparent. In fact, region 
A in NGC 3310 and region 74C in NGC 4258 show the same 
L(WR) /H/3 ratio despite their different metallicity (by a fac- 
tor of about 3). This implies that both regions have similar 
ratio of WR to O stars, while, at the same time, their effec- 
tive temperatures are considerably different. 

The interpretation of the equivalent widths of WR fea- 
tures is even more difficult, since they may be affected by 
the continuum from any underlying stellar population. Com- 
posite populations have been postulated for region A in 
NGC 3310 (Pastoriza et al 1993) and region A in NGC 7714 
(Garcia- Vargas et al. 1997). Under this assumption, the 
equivalent widths, EW(WR) and EW(Hell), of the cluster 
with WR stars are increased to 6 and 4.8 A for region A 
in NGC 3310 and 3.6 and 2.5 A for region A in NGC 7714 
which might still be compatible with Scharer & Vacca mod- 
els. In the case of region 74C, combinations of two ionizing 
clusters can be found which adequately reproduce the emis- 
sion line spectrum; however, this assumption would produce 
equivalent widths for the WR cluster higher than the maxi- 
mum predicted for the models with Z=0.008, although prob- 
ably compatible with solar metallicity models. Again, a finer 
metallicity grid for the WR cluster synthesis models would 
be needed in order to explore this possibility in further de- 
tail. 



6 SUMMARY AND CONCLUSIONS 

We have analyzed eight Hll regions in the LINER galaxy 
NGC 4258 using spectrophotometric observations between 
3700 and 9650 A. For two of the regions it has been possible 
to measure the electron temperature from the [Sill] A 6312 
A line, which allows the derivation of accurate abundances 
following standard methods. For the rest of the regions an 
empirical calibration based on the sulphur emission lines 
has been used to determine a mean oxygen content. The 
derived metallicities range from 0.3 to 0.6 Zq. In particular, 
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Table 7. Evolutionary models for region 74C 

Parameter Observations 5.3 Myr Model (GBD) 4.5 Myr Model (ST99) 
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Table 8. WR feature intensities and equivalent widths in region 74C and other extragalactic HII regions. 



Region 

Region A (NGC 3310) 
NGC 604 (M33) 
Region A (NGC 7714) 
Region 74C (NGC 4258) 



L(WR)/H/3 
0.16 
0.12 
0.07 
0.15 



EW(WR)(A) 

5.7 
10.3 

1.5 
10.8 



L(HeII)/H/3 
0.11 
0.06 
0.05 
0.09 



EW(HeII)(A) 
3.9 
5.2 
1.0 
6.6 



12+log(0/H) 
8.21 
8.40 
8.47 
8.62 



(K) 



40000 
37000 
36000 
35300 



the metallicities found for regions 74C and 5N, previously 
reported by Oey & Kennicutt (1993) to be close to solar, 
are found to be lower by a factor of two. 

For each observed region, we have also estimated the 
functional parameters: ionization parameter and the effec- 
tive temperature of the ionizing cluster. Most regions show 
ionization parameters of the order of 10 -3 and effective tem- 
peratures of around 35300 K, except regions 74C and 69C 
for which higher ionization parameters are found. We have 
also derived the physical properties of the regions and their 
corresponding ionizing clusters: fillinO factor, mass of ionized 
gas and mass of ionizing stars. Most of the regions, except 
74C, have small ionizing clusters with masses in the range 
4000 to 25000 solar masses. These values constitute in fact 
lower limits since the regions are assumed to be ionization 
bounded and the presence of dust has not been taken into 
account. 

WR features have been detected in region 74C and, to a 
lower extent, in region 5NA. A detailed modelling has been 
carried out for region 74C using different sets of models: 
those of Schaerer & Vacca (1998) for WR populations and 
those of Leitherer et al. (1999) and Garcia- Vargas, Bressan 
& Diaz (1995) for ionizing populations, to try to reproduce 
simultaneously both the WR features and the emission line 
spectrum. No consistent solution has been found: the models 
which better reproduce the WR features (a cluster 4.5 Myr 
old in the models of Schaerer & Vacc8 and Leitherer et al.) 
result too hard to reproduce the emission line spectrum. 
Conversely, a cluster 5.3 Myr old in the models of Garci'a- 
Vargas, Bressan & Di'az (1995) adequately reproduces the 
emission line spectrum, but produces too few WR stars to 
explain the observed WR features. 

When comparing the WR "bump" and Hell intensities 
and equivalent widths with other well studied HII regions 
of different metallicities, the data on region 74C are found 
to fall into the observed ranges. However, the expected de- 



crease of WR feauture intensities with decreasing metallicity 
is not clearly observed. On the other hand a decrease in the 
effective temperature of the ionizing radiation with increas- 
ing metallicity is apparent. This seems to imply that, while 
the number ratio of WR to O stars is similar for regions 
of metallicities differing by a factor of 3, the WR stars are 
cooler for the regions with higher metallicity. 

Both more observations of confirmed high metallicity 
regions and a finer metallicity grid for the evolutionary syn- 
thesis models are needed in order to understand the ionizing 
populations of HII regions. 
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Table 6. Single Star Photoionization Models for the observed HII regions 
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12 + log(0/H) 
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8.62 ± 0.11 
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8.55 


8.40 ± 0.15 
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8.65 ± 0.15 
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8.50 ± 0.15 
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8.50 ± 0.15 


12 + log(S/H) 
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6.88 ± 0.05 


6.94 


6.80 ± 0.09 


6.83 
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6.98 




6.98 




log(N/0) 


-0.98 


-1.03 ± 0.06 


-0.93 


-0.94 ± 0.08 
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log(S/0) 
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0.26 ± 0.01 


0.25 
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